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ABSTRACT 

We performed cosmological, magneto-hydrodynamical simulations to follow the evolution 
of magnetic fields in galaxy clusters, exploring the possibility that the origin of the magnetic 
seed fields are galactic outflows during the star-burst phase of galactic evolution. To do this we 
coupled a semi-analytical model for magnetized galactic winds as suggested by Bertone et al. 
( |2006 l) to our cosmological simulation. We find that the strength and structure of magnetic 
fields observed in galaxy clusters are well reproduced for a wide range of model parameters 
for the magnetized, galactic winds and do only weakly depend on the exact magnetic structure 
within the assumed galactic outflows. Although the evolution of a primordial magnetic seed 
field shows no significant differences to that of galaxy clusters fields from previous studies, we 
find that the magnetic field pollution in the diffuse medium within filaments is below the level 
predicted by scenarios with pure primordial magnetic seed field. We therefore conclude that 
magnetized galactic outflows and their subsequent evolution within the intra-cluster medium 
can fully account for the observed magnetic fields in galaxy clusters. Our findings also suggest 
that measuring cosmological magnetic fields in low-density environments such as filaments is 
much more useful than observing cluster magnetic fields to infer their possible origin. 
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1 INTRODUCTION 

Magnetic fields have been detected in galaxy clusters by radio ob- 
servations, via the Faraday rotation signal of the magnetized cluster 
atmosphere towards polarized radio sources in or behind clusters 
( [Carilli & Taylor|2002^ and from diffuse synchrotron emission of 
the cluster atmosphere (see |Govoni & Feretti 2004 ; Ferrari e t al.| 
|20081 for recent reviews). However, our understanding of their ori- 
gin is still very limited. 

At present, models for the origin of seed fields can be classi- 
fied in three main groups. In the first, a magnetic field is created in 
shocks through the "Biermann battery" effect ( Kulsru d et 31.11997] 
|Ryu et al . 1998, Mi niari et al. 2001). A subsequent turbulent dy- 
namo boosts it to the field strength observed in galaxy clusters. A 
second class of models invokes processes that took place in the 
early universe. In general, they predict that magnetic seed fields fill 
the entire volume of the universe; however the coherence length of 
the field crucially depends on the details of the models (see|Grasso 



& Rubinstein|2001|for a review). Finally, galactic winds (e.g. 



Volk 



& Atoyan 20001 or AGN ejecta (e.g. [EnBlin et al.1[T997{ [Furlan- 



etto & Loeb 2001 , and references therein) can produce magnetic 
fields and pollute the proto-cluster region. In such models, the mag- 
netic field can also originate from an early population of dwarf. 
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starburst galaxies ( [Kronberg et al.|1999) at relatively high redshift 

(2 4 - 6). 

In previous work, non radiative simulations of galaxy clus- 
ters within a cosmological environment which follow the evolu- 
tion of a primordial magnetic seed field were performed using 
Smooth-Particle-Hydrodynamics (SPH) codes ( [Dolag et al.|1999] 
'2002' '2005) as well as Adaptive Mesh Refinement (AMR) codes 
(Bruggen et al. 2005, Dubois & Teyssier||2008{ |Li et al.||2008) l. 
Although these simulations are based on different numerical tech- 
niques they show good agreement in the predicted properties of the 
magnetic fields in galaxy clusters, when the evolution of an initial 
magnetic seed field is followed. This work has also demonstrated, 
that the properties of the final magnetic field in galaxy clusters do 
not depend on the detailed structure of the assumed initial mag- 
netic field. The spatial distribution and the structure of the predicted 
magnetic field in galaxy clusters is primarily determined by the dy- 
namics of the velocity field imprinted by cluster formation ( [Dolag | 
|et al.|1999l |2002| l and compares well with measurements of Fara- 
day rotation. 

The creation of magnetic fields in shocks through the "Bier- 
mann battery effect" ( Kulsrud et al.||1997[ |Ryu et al.||1998[ l, and 
subsequent turbulent dynamo action can be followed as well as 
a prediction can be made for magnetic field values from velocity 
fields inferred in cosmological simulations ( |Ryu et ah 2008). Both 
methods predict magnetic field strengths in filaments with some- 
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what higher values (e.g. see |Sigl et al.|2004[ l than found in simula- 
tions that follow the evolution of a primordial magnetic seed field. 

Faraday rotation can be observed in several radio galaxies lo- 
cated at different radial distances with respect to the cluster center. 
Motivated by numerical simulations ( [Dolag et al.||200l) , the ob- 
served magnetic field is often modelled with a radially-declining 
field strength and a power law spectral structure. ^From such obser- 
vations, once can constrain the power law spectral index ( M urgia| 
|et al.|2004l |Govoni et al.^2006j or directly reconstr uct the power 
spectrum of the magnetic "field |vogt & Enfilin|2003|[2005) . Given 
the sparse observational data available at the moment, a degeneracy 
exists between the central value of the magnetic field and its rate of 
radial decline (see for example [Bonafede et al.|2008l|Guidetti et al.| 
|2008| ), for which detailed predictions from simulations can be use- 
ful in breaking the degeneracy. The simulations must therefore ex- 
amine different possible magnetic field origins in galaxy clusters in 
order to test the robustness of the inferred magnetic field properties. 

Recently the validity of models that produce a cluster mag- 
netic field from galactic winds has been supported by a semi- 
analytic modelling of galactic winds ( [Bertone et al.|2006l l. How- 
ever, these models are unable to predict how the magnetic fields 
produced by the ejecta of the galaxies are compressed and ampli- 
fied by the process of structure formation. Therefore, the structure 
of the final magnetic field in galaxy clusters cannot yet be predicted 
by these models. In this work we extend these studies by directly in- 
corporating the galactic outflow model in magneto-hydrodynamical 
simulations of structure formation. 

The paper is structured as follows: In section[2]we present the 
details of the cosmological setup, concentrating especially on the 
coupling of the semi-analytic model to the cosmological simula- 
tions. Details of the wind model and the seed magnetic field are 
presented in section [3T| and appendi x [A| The general results of our 
simulations are presented in section |4]and in section |5] we discuss 
our findings, particulary how galaxy clusters formed in our simula- 
tions. Finally, we present our conclusions in section[6] 



2 SIMULATION SETUP 

To investigate the evolution of a magnetic seed field produced by 
galactic outflows one needs to couple the galaxies formed in a cos- 
mological, magneto-hydrodynamical simulation with a model for 
the galactic outflows. 

2.1 Initial conditions 

We used a constrained realization of the local universe (see |Dolag| 
|et al.pOOBj and references therein), with initial conditions similar 
to those used by [Mathis et al.| p002| l in their study of structure 
formation in the Local Universe. The initial density fluctuations 
were constructed from the IRAS 1.2-Jy galaxy survey by smooth- 
ing the observed galaxy density field on a scale of 7 Mpc, evolving 
it linearly back in time, and then using it as a Gaussian constraint 
( [Hoffman & Ribak||I99I| l for an otherwise random realization of 
the ACDM cosmology. The volume constrained by the IRAS ob- 
servations covers a sphere of radius ~ 115 Mpc centered on the 
Milky Way. This region is sampled with high resolution dark matter 
particles and is embedded in a periodic box of ~ 343 Mpc length. 

^ [http : / / www ■ mpa-garching . mpg . de/galf orm/[ 
data_vis/ index. shtml#moviel2 




Figure 1. Shown is a visualization of the constrained local universe simula- 
tion at redshift z = using SPLOTCH jPolag et al.|2008} . Structures which 
can be identified with their counterparts in the real universe are fabefled. An 
animation flying through the simufation can be downioaded from the MPA 
WebsiteE] 

The gravitational softening length used is e = 14 kpc, which corre- 
sponds to the smallest SPH smoothing length reached in the dense 
centers of halos. In low density regions, the resolution is lower. The 
region outside the constrained volume is filled with dark matter par- 
ticles at fower resolution, ailowing a good coverage of long range 
gravitational tidal forces. Many of the most prominent clusters ob- 
served locally can therefore be identified directly with halos in the 
simulation, and their positions and masses agree well with their 
simulated counterparts. Figure [T] shows a rendering of the simula- 
tion at redshift z = Q where some of the identified structures are 
labelled. 

In this work we use extended initial conditions where the 
original high resolution dark matter particles are split into gas 
and dark matter particles with masses of 0.69 X 10® Mq and 
4.4 X 10® M© respectiveiy. The most massive clusters in our simu- 
lations are hence resolved by nearly one million particles. The grav- 
itational force resolution (i.e. the gravitational softening length) of 
the simulations was set to be 14 kpc, which is comparable to the 
inter-particle separation of the SPH particles in the dense centers 
of our simulated galaxy clusters. 

2.2 Seeding strategy 

In principal one would like to follow the dynamics of galactic out- 
flows driven by stellar activity (either star-formation or star-burst) 
self-consistently in cosmological simulations. However, so far this 
it not possible. Although cooling and star-formation processes can 
be followed in standard cosmological simulations, the detailed in- 
teraction with magnetic fields leads to regions in which the mag- 
netic field pressure exceeds the thermaf one, especially in high res- 
olution studies like those we pfan to perform. In absence of any dis- 
sipative process when performing ideal MHD simulations, these re- 
gions then dominate the time-stepping in the simulations and cause 
them to stall. We therefore took a simpler approach by coupling 
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Figure 2. Newly formed haloes per Mpc cubed as a function of redshift in 
the reference simulation, involving a primordial magnetic seed I Dolag et aL 
|2005^ . The black curve shows the cumulative distribution, the red curve the 
differential one (per Gyr). The dashed line marks the time selected for the 
single seeding, the dotted line marks the additional seeding times in the 
multi-seed run. 



a semi-analytical recipe for star-burst-driven, magnetized outflows 
( [Bertone et al.|2006) with non-radiative, cosmological MHD simu- 
lations. In general, this should be implemented as a continuous pro- 
cess, starting at very high redshift z ~ 10 ^Kronberg et al.|1999| l. 
For practical reasons, we approximate this process by an instanta- 
neous magnetization of all galaxies inferred in the simulations at a 
selected instance in time (e.g. z — 4.1), and then evolve the simula- 
tion with magnetic seed fields until the present day. As the cosmo- 
logical simulation evolves, the magnetized gas is stripped from its 
host galaxies as it falls towards the dense gas contained in galaxy 
clusters. The intra cluster medium (ICM) is thereby enriched by 
the stripped, magnetized gas of several hundreds of galaxies. Fur- 
thermore, the magnetized gas is processed by shear-flows and tur- 
bulence within the cluster atmosphere, which eventually leads to 
the magnetic field strength and structure observed today in galaxy 
clusters. To justify our approximation we also performed one sim- 
ulation where we repeated the seeding process several times and 
demonstrated that for cluster magnetic fields these further seed- 
ings are irrelevant, as they mainly involve newly-formed galax- 
ies outside the proto-cluster region. The redshift z = 4.1 of the 
single seeding event was chosen as a compromise between com- 
plete formation of all relevant galaxies (e.g. choosing a low seeding 
redshift) and avoidance of galaxy mergers or destruction in proto- 
clusters (e.g. choosing a high seeding redshift). Figure [2] shows the 
number density of formed halos as a function of redshift (differ- 
ential per Gyr and cumulative) as inferred from our reference run 
with a cosmological seed field ( [Dolag et al.|2005 V The dashed line 
marks the preferred redshift of seeding which is chosen to ensure 
that most of the galaxies in the proto-cluster region are formed and 
not yet merged with others, destroyed or stripped from their hot 
gas haloes in the proto-cluster. The dotted line marks the additional 
seeding times used in the multi-seed reference run. As seeding tar- 
gets we choose galaxies which are newly formed (i.e. since the last 
seeding) and whose mass is smaller than A/naio = 3 x 10^^ Mq 
since our wind model does not apply for group like objects. Note 
that the constraint on mass is only of relevance for the low redshift 
seeding in the multi-seed reference run. 



tsb [yr] 

Figure 3. Evolution of the magnetic energy in a wind bubble of a small 
galaxy (Mism = 10 x W^^Mq) for different disc and halo parameters 
(Bq, Bq). The dashed line marks the star-burst time of 150 Myr jde Grijs| 
|200H , the dot-dashed lines mark the envelopes of the two wind phases. For 
the chosen t^^, the initial energy in the wind bubble Bq does not influence 
the bubble energy at late times. 

2.3 Numerical Method 

All simulations were performed using the MHD implementation in 
GADGET-2 ( |Springel||2005 t, similar to that used in 'Dolag et al. 
1 2005 ^. Detai ls of the MHD implementation can be found in Dolag 
|& Stasyszyn|( |200^ . The code is based on an entropy-conserving 
formulation of Smooth Particle Hydrodynamics (SPH) l |Springel"&| 
|Hemquist|2002) . It was supplemented with a treatment of magnetic 
field using ideal MHD (see |Dolag & Stasyszyn|2008| >. Besides fol- 
lowing the induction equation for the magnetic field, we take mag- 
netic back-reaction into account using a symmetric formulation of 
the Lorentz force based on the Maxwell tensor. The treatment of 
magnetic fields in SPH was improved by explicitly subtracting the 
part of the magnetic force which is proportional to the divergence of 
the magnetic field, as described in B0rve et al. ( 2001 ). This helps to 
keep numerically induced divergence of the magnetic field at negli- 
gible values. It also helps to avoid instabilities of the MHD formu- 
lations in SPH, especially in regions where magnetic field pressure 
substantially exceeds the thermal pressure. 

Note that in simulations without radiative cooling like ours 
the magnetic field pressure stays well below the thermal one, even 
within the cores of the most massive galaxy clusters. In very strong 
shocks however it can still happen that the magnetic field is com- 
pressed so substantially that magnetic forces dominate the thermal 
ones for brief periods of time. Such situations are handled more 
accurately with our new formulation, which includes several other, 
small improvements in the numerical treatment. For detailed in- 
vestigations and information on the code and its performance see 
[Dolag & StasyszynlpOOS) . 



3 SIMULATIONS 

To identify galaxies in our simulations, we applied SUBFIND 
( [Springel et"d:i|200T) which uses a friends-of-friends (FoF) 
algorithm to identify locally over-dense, self-bound particle groups 
associated with galaxies, even when they are inherited within a 
larger parent group. Originally, this algorithm was based on pure 
dark matter simulation, so we utilized a modified version to handle 
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Parameter 



Value 



Source 



Bo 
Bg 

tsh 

MiSM 



400 pc 
50, 5, 0.5 
30, 3, 0.3 

lOMe/yr 
150 Myr 

< 3 X lO^^Mp- 



|Kleinetal.|l988| 
(Klein etal. 19881 
(Chyzy & Beck 2004] , 
^Soida et al. 20Ul| 
(de Grijs 2001 
(de Grijs 20UT 
from simulation 



Table 1. Summary of the parameters used for the wind model. Except for 
Bq all other parameters are based on observations of M82. The parameter 
set is scaled to the simulated galaxies by Mism 



the gas component as well (for more details see Dolag et al. 2008). 



3.1 The Adopted Wind Model 

To obtain the size and strength of the magnetic field in a star- 
burst driven galactic outflow for each of the identified galaxies, 
we adapted the semi-analytical model by |Bertone et al.| p006|l. 
This model assumes an adiabatic expansion of a spherical gas bub- 
ble with homogeneous energy density, fed continuously by the 
starburst-driven outflows. The initial bubble size (before the star- 
burst) can then be characterized by a field strength Bo and a radius 
Rq. The bubble expansion is driven by the starburst in the galaxy, 
which is expelling gas at the rate Af*, dragging the frozen-in mag- 
netic field Bq into the ICM. The wind velocity is a function of the 
star formation rate, following |Shu et ar| ( |2005| l we set the ISM con- 
stant K = 0.5, which determines the dependence of wind velocity 
and mass outflow rate on ISM properties (evaporation parameter 
and blast wave speed): 



/ . \ 0.145 

The galaxy thus injects magnetic energy into the wind at 



(1) 



a rate 



' Pin 



where the injected mass rate is ( |Shu et al.|2005| l 

/ \ 0.71 



yr 



(2) 



(3) 



and the average injected mass density pin follows from the injected 
mass blown through a spherical surface at the galactic radius, which 
is assumed to be a fraction of the virial radius Rg — 0.1 -R200- 



Pin 



(4) 



The injected magnetic field ener gy density eg^^ derea ses adiabati 
cally as B oc pa (see App.A m 



Bertone et al. 



20061): 



(5) 



where pisM = Mjsm/ ^Rq is the average ISM mass (Mism) den- 
sity inside a sphere of galactic radius. Finally, the predicted injected 
magnetic energy will be 



Neglecting shear amplification the time evolution of the energy in 
the sphere is given by 



dt 



1 V^t) 
3 l/w(t) 



Esit) 



(7) 



(8) 



Here local alignment between magnetic field and wind direction 
can not be assumecffjbecause of turbulent instabilities and random 
motions in the gas. Still global structures might develop, because 
the magnetic field is dynamically unimportant and follows the pre- 
ferred wind direction, parallel to the rotation axis of the halo. Fig- 
ure [3] shows the evolution of the bubble energy as a function of 
the effective starburst time for several different values of the disc 
and halo parameters (Bg, Bq) of the model. The model is based on 
a galaxy with ISM mass of AfisM = 10 x IQ^'^Mq. The initial 
bubble energy is diluted on reasonable starburst time scales, and 
therefore does not affect the final magnetic energy contained in the 
wind. The starburst time t^h is therefore degenerate with Bq over a 
wide energy range. For more details on the wind model see |Bertonel 
[etaL] ( |2006) . 

Given the mass-cut on A/ism, the resulting wind velocity is suf- 
ficient to reach every particle of a halo during the star-burst time 
scale assumed. Therefore we always seed the complete halo, by 
setting the magnetic field for every particle. Due to the friends- 
of-friends algorithm used to identify the structures, overlapping of 
halos is impossible and a minimum halo distance is defined by the 
linking length. 



3.2 Applying the Wind Model 



Contrary to |Bertone et al.| ( [2006^ , we cannot follow (and integrate) 
the evolution of the wind model over cosmic time based on semi- 
analytical modeling. Instead we are holding the numerical simula- 
tion at certain epochs and identify newly formed galaxies. We then 
integrate the wind model for every of these galaxies assuming a 
generic star-formation rate and a generic star-burst time. This ef- 
fectively mimics the magnetic seeds obtained from such a galaxy 
at these epochs. The procedure is applied (in an approximate way) 
instantaneously. Given the mass-cut on A^ism, the resulting wind 
velocity is so large that the wind can reach every particle within 
the virial radius of the halos during the assumed star-burst time 
scale. Hence, we are always seeding all particles within the virial 
radius, which by construction do not overlap. Continuing the mag- 
netohydrodynamic simulation we follow the complete dynamics of 
the magnetized atmospheres, including their stripping and mixing 
in the denser environment (e.g. filaments or clusters atmospheres). 
The galactic virial radius selected for the seeding procedure is typi- 
cally smaller than the radius of 100 kpc inferred for the wind driven 
bubbles in [Bertone et al.| ( [2006[ l. The magnetised bubbles of the 
galaxy sized halos are stripped by ram pressure effects, which are 
already important in filaments (see ( [Dolag e t al. 2006), ( [Saro et aT] 
[2006[ l, ( [Dolag et al.[[2008] l), as well as partially ejected by inter- 
actions with other galaxies. The stripped gas gets mixed with the 
cluster atmosphere, and large scale motions typically distribute this 
magnetized material across a much scale than the expanded bubble 
size predicted by the semi-analytical models. 

The instantaneous seeding process is implemented in the sim- 
ulations in the following way: 



^Bi„ 



\BlRlv. 



(6) 



This results in the factor 1/3 in the dilution term 
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• The simulation is stopped at the seeding redshift (usually z = 
4.2), and the newly formed galaxies (since the last seeding) are 
identified. 

• The gas mass of the halos identified in the simulation is com- 
bined with the wind model to estimate the magnetic energy con- 
tained in the wind (see section [3T| l. 

• The magnetic moment of a dipole/quadrupole having the same 
magnetic energy is calculated. A smoothing length of 14kpc is 
used as softening length in the integral (appendix [Ajl. 

• The magnetic field vector of every gas particle in the corre- 
sponding halo is set to match the dipole/quadrupole. The field of 
every particle inside the softening length is set to a random orien- 
tation and the field strength at the softening length. 

• The simulation is continued with the newly added magnetic 
field. 

The parameters of the model can in principle be inferred from 
observations of individual galaxies. Table [T] shows some observa- 
tional constraints on the parameters as obtained from one of the 
best observed, starbursting galaxies, M82. Note that in starbursting 
galaxies usually disc fields up to 100 /iG are observed (R.Beck, 
priv.com.). This is the case on scales smaller than the resolution 
archieved by our simulations (14 kpc). We rather follow the con- 
servative approach in |Bertone et al.| ( ["2006^ using a field of several 

tens of ^G. 

M82's observed mass of 6 X 10^ Mq ISofue' 1998) is smaller than 
the average mass of the seeded halos (table [2](. To take this into 
account, we assume that the mass of the ISM Mism in the model 
is equal to that of the halo mass in the simulation, and not from 
M82. There is further mass dependence in the mass outflow rate 
M*. Therefore the value inferred from M82 might be underesti- 
mated and the obtained wind energy represents a lower limit for 
large halos. 

As the values quoted are still quite uncertain, we performed sim- 
ulations where we varied the value for the halo and disc magnetic 
field (_Bo, Be) within a wide range. 

The magnetic field structure within galactic outflows and 
the properties of the resulting wind driven bubble are largely 
unconstrained as well. In M82, |Reuter et al.| ( |1994[ l observe a 
symmetric poloidal magnetic field structure in the core, while in 
edge-on galaxies combinations of symmetric (SO) and antisym- 
metric (AO) field structures are observed (R.Beck, priv.com.). AO 
and SO dynamo modes can be approximated by a softened dipole 
and quadrupole field structure, respectively. We therefore describe 
the structure of the magnetic seed field as a (softened) dipole field, 
which we normalized so that the field energy corresponds to the 
energy in the magnetic bubble inferred from our semi-analytical 
model. The orientation of the dipole is chosen to align with the 
spin of the underlying dark matter halo. To verify the dependence 
of our results on the detailed structure of the seed fields we also 
performed a run where we used a quadrupole instead of a dipole 
structure for the magnetic seed field. 

Including the reference run (Control Run) which follows the 
evolution of a cosmological seed field, we performed six different 
simulations to explore the parameter space. In all runs the starburst 
time tsb is kept fixed at 150 Myr. We performed three simulations 
for different values of the magnetic fields (halo Bo and disc Bq ) 
which all have a dipole like structure. Starting from the M82 like 
value of 50/iG {Dipole), we also performed a run with one tenth 
of the value (0.7 Dipole) and one hundredth {0.01 Dipole). We 
used the Dipole parameter set from above to perform a run with 
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Table 2. Halo statistics of the multiple seed run. The rows give the output 
number, the corresponding redshift, the number of new haloes identified 
since last seeding and the mean mass of these haloes, respectively. 

a quadrupole structure {Quadrupole^. We further repeated the 0.1 
Dipole run and performed four additional seeding episodes {Multi 
Seed) to all newly formed haloes as summarized in table [2] Form 
the results (e.g. compare left and right panel in the middle row of 
figure |5j its clear that after the first seeding, many new haloes are 
still seeded, but they are typically smaller in mass (see last column 
in table [2]( and preferentially located outside the proto-cluster re- 
gion. 



4 GENERAL RESULTS 

Figure |4] shows the results from the two different seeding strate- 
gies compared in this work. The right column shows the case of 
primordial magnetic seed field {Control Run) and the left column a 
run which used the seed fields of galactic outflows {Dipole). Shown 
is the magnetic field amplitude projected through the box. The top 
row shows the field at the first seeding of the galactic outflows while 
the lower row shows the final magnetic field at redshift zero. One 
can clearly see the differences between the two models at early 
times: the low-level primordial magnetic field fills the whole vol- 
ume, whereas the magnetic field from the galactic outflows is con- 
centrated on the galaxy population. At low redshift the magnetic 
field in galaxy clusters appears comparable in both models. This 
reflects the fact that inside the galaxy clusters the magnetic field 
is strongly processed by compression, shear flows and turbulence. 
Therefore, the final magnetic field is shaped by such processes 
rather than by the initial conditions, in agreement with previous 
findings jPolag et al.|1999||2002^ and analytical mo deling of satu- 
rated dynamos in clusters ( |Subramanian et al.|2006[|EnBIin & Vogt| 
|2006[ l. However, the magnetic field in filaments looks quite differ- 
ent. In the model with galactic outflows the magnetic field in the 
collapsed objects along the filaments is higher than for a primor- 
dial seed field; the magnetic field inside the diffuse component of 
the filaments is much more prominent in the case of cosmological 
seed fields, 

because the density in the filaments is too low to strip the hot 
and magnetized atmosphere of the galactic halos completely. Note 
that our models do not include a kinematic component which could 
lead to an evaporation of the haloes around galaxies. This and espe- 
cially the effect of AGN-driven outflows could change the picture, 
as the outflows could fill volumes of Mpc size in these low den- 
sity environments (see |Kronberg| ( (2006) , and references therein). 
We also note that in principle small galactic halos (down to masses 
of 10*A/q) - which are not resolved in the simulation - could 
contribute to the magnetization of the low density environment. 
However, observations suggest that the star-formation rate of such 
galaxies would be quite small (far below the value of lOM0/yr 
Dad di et al.| |2007| adopted in this work for resolved galaxies). 
Wence, we do expect them to contribute significantly. 
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Figure 4. Visualisation of tlie magnetic field strength in the simulation box at redshifts z = A.l and z = 0. Every image shows a region of a linear size of 
204 Mpc/(1 + z) and was made using the same color bar. The upper left panel shows the magnetic field due to instantaneous seeding with a dipole M82-like 
structure at this redshift with maximum field strengths of Ri 5 nG. The upper right shows a simulation with homogenous cosmological magnetic seed field 
jPolag et al.|2005f . Here, the field strength more continuously traces the underlying structures of the matter distribution and reaches values of up to ^ 10 nG 
in the highest density regions. The lower two panels show both simulations (left Dipole, right Control Run) evolved to z = 0. 



Figure[5]shows the magnetic field amplitude projected through 
the box for all different runs at redshift zero. When the halo mag- 
netic field in the wind model is reduced, the field within the indi- 
vidual galaxies is strongly reduced. Therefore the magnetic fields 
within the filaments become less prominent (Dipole, 0.1 Dipole 
and 0.01 Dipole panels). However, the magnetic field within the 
collapsed structures of galaxy clusters only slightly changes, indi- 
cating that the amplification within the cluster centers is strongly 
influenced by saturation effects, well in agreement with previous 
findings ( |Dolag et al.||2005) . The situation changes for the multi 



seed run {Multi Seed). Although again the magnetic field in the 
clusters is not affected, a clear change is visible in the filaments, 
where many more seed fields from galaxies appear. This is be- 
cause galaxies are formed earlier in high density environments like 
proto-cluster regions than in low density environments, like fila- 
ments. Subsequent seeding of newly formed haloes therefore pref- 
erentially happens in low density environments like filaments. 

A more quantitative way of displaying the differences of our 
models is to show the volume weighted cumulative filling factors 
obtained within the whole box. In figure|6]the rising (falling) curve 
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Figure 5. Visualization of the magnetic field strength in the simulation box at redshift 2 = 0. Every image shows a region of 204 Mpc/(1 + z), using the 
same arbitrary color bar. Shown are the results of the Dipole (top left), 0.1 Dipole (middle left), O.OI Dipole (bottom left), Quadrupole (top right), Multi Seed 
(middle right), and the Control simulation (bottom right), respectively. 
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Figure 6. Volume-weighted cumulative filling factor of the magnetic field 
in the cosmological box. The filling factor is plotted above and below a 
threshold as a function of the threshold (the rising and falling curves). The 
black lines show the cosmological seed field from Dolag et al. 1 2005}. Also 
shown are the wind-induced seeding from the Dipole (red) and Quadmpole 
(blue) and the curves from the simulations with decreased field strength 
(green curve: 0.1 Dipole; magenta curve: 0.01 Dipole) . The results from 
the simulations with the Multi Seeds are given by the light blue curve. 



indicates the relative filling of the box by a magnetic field weaker 
(stronger) than a certain value, for all simulations at redshift z = Q. 
In the very high field regime (/iG and above) all simulations behave 
similiar to each other, again indicating that amplification within 
galaxy clusters is strongly influenced by saturation effects. In the 
intermediate regime (most visible in the rising curves), the simu- 
lations split in three groups, distinguished by the chosen value for 
the disc magnetic field, which controls the amount of magnetic field 
within the galactic outflow. In this regime, magnetic field amplifi- 
cation reflects the underlying velocity field structure. The simula- 
tion, which evolves a cosmological seed field {Control Run), cor- 
responds (by chance) to a medium value of the galactic fields {0.1 
Dipole). At low magnetic field (best visible in the rising curves), 
the filling factor is dominated by the properties of the seed field. 
Therefore, all models are clearly distinguishable. This region cor- 
responds to the low density regions, including the diffuse matter in 
filaments. For lowest fields, the four single seeded simulations con- 
verge to a non zero value, representing the un-magnetized volume 
left in the box. The multi seeded simulation shows a lower fraction 
of unmagnetized volume, which is expected from the seeding strat- 
egy. In the simulation using a cosmological seed field, no volume is 
filled with a field smaller than this seed field (adiabatically diluted 
with the local density). 

This shows that for galaxy clusters the single seed simplifica- 
tion in our approach is not crucial. On the other hand, low density 
regions like filaments provide a natural testbed for the origin of 
cosmological magnetic fields and also are very sensitive to the ap- 
proximations and simplifications made. 

This is expected from previous findings. In simulations where 
a cosmological magnetic field is evolved, the field amplification is 
linear in low density regions undergoing mainly adiabatic compres- 
sion or dilution, whereas the orientation is typically aUgned with 
the filament ( Dolag et al.| pOOSl i; [Briiggen et aL| pOOSl l). There- 
fore, such regions are only mildly influencing the magnetic field 
configuration imposed by the seeding mechanism. 



5 CLUSTER BASED ANALYSIS 

The amplification of magnetic fields in galaxy clusters is driven 
by adiabatic compression, and merger events, which induce shear 
flows and subsequent turbulent flows. All these processes depend 
on the cluster environment and therefore may depend on cluster 
mass. Observationally there is yet not enough data available to 
track how the magnetic field varies within clusters as function of 
radius or how the mean field within the cluster varies with cluster 
mass (or equivalent, with the mean cluster temperature). Cassario] 
|et al.||2007| concluded from the properties of observed, diffuse ra- 
dio emission (so-called radio halos) that the mean magnetic field 
within the radio emitting region appeared independent of cluster 
mass. 

But the volume of the radio-emitting region is found to in- 
crease with cluster mass. 

Therefore, even if the central magnetic field does increase with 
cluster mass, a constant mean magnetic field could still be inferred, 
if the shape of the radial field profile decreases properly. Further- 
more, although indications for a radial dependence of the magnetic 
field can be inferred from clusters where rotation measures of sev- 
eral radio galaxies are available ( [Dolag et al.|200l] [Murgia et al.| 
|2004[ |Govoni|2006^ , there exits still a large degeneracy between 
the parameters describing the magnetic field in clusters ( Bo nafede| 
|et al. 2008; Guidetti et al.i2008). Hence, it is quite important to 
compare observed scalings with predictions from simulations. Es- 
pecially the dependence of scaling relations on the adopted model 
for the origin of the magnetic seed fields has to be studied in more 
detail. To this end, various simulations performed in the past gave 
quite different results. 

Models based on shock induced field generation predict mag- 
netic field strengths with a very weak dependence on the distance 
to the cluster center a nd also predic t a scaling of the mean magnetic 



field B oc VT (Miniati et al. 



2001 



In contrast, previous SPH simulations following a primordial 
magnetic seed field predict a steep radial dependence for the mean 
magnetic field, mainly following the decline of the density ( |Dolag| 
|et al.| 200r, 2005 1. Such steep profiles were confirmed with AMR 
simulations ( Briiggen et al. 2005, Dubois & Teyssier_2008) . SPH 
simulations also predict a very steep correlation of th e mean mag 



netic field with the cluster temperature, e.g. B oc T i Dolag et al. 
,1999^ 



5.1 Radial profiles 

In figure |7] we plot angularly averaged radial profiles of density 
(left column) and mass-weighted temperature (right column). The 
individual lines are taken from the 16 most massive clusters within 
the simulation. They are scaled to the virial radius of the individual 
clusters and are normalized to the same mean value within 0.1i?vir- 
The upper row shows the results from the control run following a 
primordial magnetic seed field jDolag et al.|2005 \ the lower row 
shows results from the dipole run with the largest value for the 
galactic halo field. Although individual haloes show small varia- 
tions because of the slight changes magnetic fields introduce in the 
systems dynamics (by adding small perturbations of the force field 
at high redshift), there is no significant dynamical influence of the 
magnetic field. This is in good agreement with the relatively small 
values of the magnetic fields (< lO^G) within the simulated galaxy 
clusters, and also with previous results. 

Figure [8] shows the RMS magnetic field profiles for the 16 
most massive clusters for all our simulations. As before, they are 
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Figure 7. Density (left column) and mass weighted temperature (right column) profiles for the 16 most massive clusters in the simulations. Here, only the 
Control Run and the Dipole Run are shown, demonstrating the absence of a significant dynamic influence of the magnetic field in the models. 



scaled to the virial radius of the individual clusters and are normal- 
ized to the same mean value within 0.1i?vir- In good agreement 
with previous findings, the magnetic field profiles generally follow 
a very similar decline as the density, but with somewhat more scat- 
ter. However, the scatter of the individual magnetic field profiles 
seems to depend on the overall strength of the magnetic field. The 
simulations with higher magnetic seed fields (e.g. the Dipole and 
Quadrupole runs) show significantly less scatter. This is a clear in- 
dication that in these simulations saturation effects dominate the 
amplification of the magnetic field in outer parts. Therefore, the to- 
tal field amplification is less dependent on the system's dynamics 
and so the differences in the scaled profiles are reduced. There is 
also some indication that the scatter in the profiles is slightly in- 
creased for the multiple seed run, indicating that there could be still 
some influence of the late forming galaxies on the outer parts of 
the magnetic profiles. In general, aside from the scatter in the pro- 
files the exact details of the origin of the magnetic seed fields do no 
strongly influence the predicted shape of the radial profiles. There- 
fore, the shape of the magnetic profiles is a robust prediction of 
the simulations, driven mainly by the dynamics imprinted from the 
structure formation process. 



dynamical state is very important. Here, the magnetic field amplifi- 
cation is still in a regime, where individual merger events strongly 
contribute to the amplification of the magnetic fields. At intermedi- 
ate and high temperatures the scatter reduces and the trend becomes 
very regular. For strong magnetic field models the mean field flat- 
tens horizontally as a function of temperature, which indicates a 
saturation of the magnetic field amplification. This can be seen well 
when comparing runs with varying dipole strength (left column of 
figure |9j. Again, there is no significant dependence on the origi- 
nal magnetic field structure assumed within the galactic outflows 
as evident from comparing the 0.1 Dipole with the Quadrupole 
simulation. Multi Seed simulations only affect low mass systems. 

The normalisation of the average magnetic field strength 
might be sensitive to the resolution of the simulation, as well as to 
the exact details of saturation of the amplification mechanism. This 
includes the final amplitude of the mean magnetic fields as well 
as the temperature (mass) scale at which clusters magnetic fields 
reach the saturation regime. The central magnetic fields (several 
^G) found in our simulated clusters are in good agreement with 
observations (e.g. see val ues ofjGuidetti et al.|2008||Govoni et al.| 



|2006|putgia et al.|2004[|Bonafede et al.|2008| > 



5.2 Average field strength 

Figure |9] shows the mean magnetic field within one tenth of the 
virial radius as a function of mean, mass-weighted temperature for 
the 16 most massive clusters from our simulations. There is a strong 
trend of the magnetic field with temperature (e.g. mass). At low 
temperatures (e.g. low masses) a large scatter is present, not only 
for models with different seed fields, but also for clusters with simi- 
lar temperature. Simulations of low mass/temperature systems also 
show a quite steep relation, indicating that for small systems, the 



5.3 Synthetic rotation measurement profiles 

To compare the simulated magnetic field from different seeds with 
observations, we repeat the comparison with observational data 
as m Polag et al.|[2002l |2005}. Fi gure [To] shows a combination 
of three observational samples, measuring the Faraday rotation of 
point like sources in or behind Abell clusters as function of their 
distance to the center of the clusters jKim et al. |1 99 1| [Clarke et aT] 
[200T][Johnston-Hollitt & Ekers|2004| >. The black line shows the re- 
sult obtained for the median of the absolute values when radially 
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Figure 8. Radial profiles of the RMS magnetic field strength for the 16 most massive clusters. The profiles are scaled to i?vir and normalized to the same 
mean value within 0.1 X i?vir- Shown are the results of the Control simulation (bottom right), Dipole (top left),0.i Dipole (middleleft), 0.01 Dipole (bottom 
\eft),Quadrupole (top right) and Multi Seed (middle right), respectively. 



binning using bins of 15 data points each. All three individual sam- 
ples are statistically compatible with each other and we combine 
all of them to allow a finer binning of the data. To calculate error 
bars we used the RMS of the median obtained by bootstrapping 
the samples of each bin 1000 times. We over plotted the values 
inferred from three elongated sources (triangles) observed in the 
single galaxy cluster A119 (Feretti et al. 19991 and one elongated 
source within the Coma cluster (diamonds) (Feret ti et al.||1995l ). 
As discussed in previous work jPolag et"al.|2005| l, due to the con- 
struction of these observational samples, the underlying selection 
function for the contributing galaxy clusters is ill-defined. This is 
especially important, as the signal depends on the mass of the un- 



derlying galaxy clusters, both due to the larger line of sight contri- 
bution from massive systems as well as the predicted dependence 
of the magnetic field strength on the cluster mass. One therefore 
expects the median profile to depend on the mass function of the 
selected clusters. However, as the cluster sample is composed of 
Abell clusters, it is reasonable to assume that it mainly consists of 
comparatively massive galaxy clusters. The additional data points 
inferred from the two massive systems Al 19 and Coma (black sym- 
bols) follow the curve from the combined sample (black line) res- 
onable well, indicating that the observed sample indeed compares 
well to the predictions for more massive systems. 

From the simulations we calculated synthetic Faraday rotation 
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Figure 9. Same amngement as figure[8]showing the mean magnetic field strength inside one tenth of the virial radius as function of the mean, mass-weighted 
temperature within the same radius. Clusters that have an observed counterpart are marked with different colors. 



maps for the 16 most massive clusters. Binning the individual maps 
in radial bins we computed the median of the absolute value of the 
Faraday rotation combining the same radial bins from all clusters. 
This procedure was repeated for two subsets of the clusters, where 
we restricted the sample to clusters with masses above 3 x IO^^A/q 
and 5 x 10^'' Af© respectively. The results are shown as three col- 
ored lines in figure [TO] In general, the shape of the resulting radial 
profiles compares to observations quite well; there is in particular 
no noticeable difference between the Control run, which follows 
the evolution of a primordial magnetic field and the Dipole 0.1 run. 
The Dipole 0.1 fits also best regarding the different strength for the 
galactic halo field. In agreement with results presented before, there 
is no visible difference between the quadrupole and the dipole con- 
figurations for the magnetic field in the galactic winds. There are 



also only very mild effects visible in the outer parts of the Faraday 
rotation profiles for our run with the multiple seeding episode. 

Finally, we calculated the projected structure function from 
some of the synthetic Faraday rotation maps. This was done follow- 
ing |Murgia£tar]([2004|, who calculated this for observed Faraday 
rotation maps: 

S'^'^\dx,dy) ={\RM{x,y)~RM{x + dx,y + dy)\), (9) 

with dx and dy being the offsets from a pixel at position {x,y). 
The resulting matrix is then averaged in radial bins to finally ob- 
tain the structure function. In figure[TT|we show the structure func- 
tions obtained from three of our clusters normalized at the largest 
scales. For the different models of the magnetic seed field the struc- 
ture functions only reveals some small differences in the Faraday 
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Figure 10. Same arrangement as figure[8]but showing the radial profiles of Faraday rotation. We binned the absolute value of the observed Faraday rotation 
as function of distance to the centres of the clusters and plotted the median forcing the radial intervals to contain always 15 data points. The data are obtained 
from combining three samples based on Abell clusters jKim et al.|1991[[Clarke et al.|200ll|Johnston-Hollitt & Ekers|2004^ . The error bars are obtained by 
bootstrapping the data points within each bin. We also included in the plot the values inferred from three elongated sources (triangles) observed in the single 
galaxy cluster A119 ^Feretti et al.|1999) and one elongated source within the Coma cluster (diamond) jFeretti et al.|1995) . For the simulations we build the 
median of the Faraday rotation for the 16 most massive clusters (red line) and subsets restricting their mass to be larger thanS X and 5 X W^'^Mq, 

blue and green line respectively. 



rotation of the first cluster. In general, there are no significant in- 
dications of the seed magnetic field in the final Faraday rotation 
structure, especially given that the effective resolution limit of the 
simulation is about 15 — 20 kpc in the central parts of the clusters. 



6 CONCLUSIONS 

In this work we performed cosmological, magnetohydrodynamic 
simulations following the evolution of magnetic fields on large 
scale structures and in galaxy clusters. Coupling a semi-analytic 
model for magnetized galactic winds as suggested by |Bertone et al.| 
((2006) with our cosmological simulation we explored the possibil- 
ity that the magnetic fields in galaxy clusters originate from galac- 
tic outflows during star-burst phases, further processed by struc- 
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Figure 11. Shown is the structure function (for details see text) calculated 
from synthetic Faraday rotation maps obtained from different clusters dis- 
played in the individual panels. The different lines correspond to different 
magnetic seed fields. 



• The general shape of the predicted Faraday rotation profile 
within clusters compares well with the sparse observational data 
available. Models that assume a field strength of 5/iG within the 
galactic halo reproduce the observed Faraday rotation profiles bet- 
ter than models with a ten times stronger or a ten times weaker halo 
magnetic field. 

• The properties of the final magnetic field in galaxy clusters 
do not depend on the exact field configuration within the magnetic 
outflows. This confirms previous studies that the structure of the 
magnetic field in galaxy clusters is primarily driven by the velocity 
field induced by the structure formation process. 

• In massive galaxy clusters, the magnetic field amplification 
saturates around values of several fiG. The mass (or temperature) 
scale on which this happens depends on the strength of the mag- 
netic seed field, and probably also on the resolution of the simula- 
tion. 

• In systems where saturation effects start to play a significant 
role, we observe only a small scatter in the magnetic scaling rela- 
tions, and in the shape of the radial, magnetic profiles. 

• In clusters where saturation effects are negligible the strength 
and configuration of the magnetic field strongly depends on the dy- 
namic state. Therefore, we observe a large scatter in the magnetic 
scaling relations and in the shape of the radial, magnetic profiles. 

• Within galaxy clusters, the structures predicted from synthetic 
Faraday rotation maps do not depend significantly on the magnetic 
seed field, and agree well with observed ones. 

• In low density environments imprints of the magnetic seed 
fields are still present and can be observed, in principle. In these en- 
vironments the galactic outflows forming at later epochs contribute 
significantly to the magnetic field configuration. 

In summary, we confirm that galactic outflows and the subse- 
quent action of structure formation can explain the properties of the 
magnetic fields observed in massive galaxy clusters. Their strength 
and shape is attributable to the velocity field induced by structure 
formation. We find that there are no measurable imprints of the 
magnetic seed fields left in the synthetic Faraday rotation maps of 
our simulated clusters. However, low density environments, like fil- 
aments, still contain significant information on the magnetic seed 
fields. Therefore, they may be used to discriminate proposed sce- 
narios once information on magnetic fields within these regions be- 
comes available with the next generation of radio instruments. 
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ture formation. We compared our results with the ones obtained 
by following a primordial magnetic seed field pdag et al.|2005l l. 
Performing several simulations, we explored the effect of vari- 
ous parameters of the adapted semi-analytic model relevant for the 
strength of the magnetic seed field from the galactic outflows. We 
also explored the effect of the magnetic field configuration assumed 
for the galactic outflows and of the seeding strategy. Our general 
findings are: 

• The typical magnetic field strengths of several fiG in galaxy 
clusters as obtained from observations of Faraday rotation are well 
reproduced for a wide range of parameters of the galactic outflow 
model. 



APPENDIX A: CONVERSION FORMULAE 
Al Dipole Energy 



E- 



B^dV 



B 



(Al) 
(A2) 



1 

1 3n (n ■ m) ~ m 

47r \r\ 

To calculate the magnetic moment m of a dipole of magnetic 
energy Eb, one assumes, that m| |ez, as the moment can be aligned 
into any other direction by a simple rotation of the coordinate sys- 
tem. The unit vector ft then becomes zj \ z\. 
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Since the field B diverges for r ^ we use a softening lengtli 
of e = 14, which corresponds to the SPH smoothing length. The 
outer boundary -Rmax takes the finite size of the halo into account. 
Setting Ho = lit follows: 

_ ^ f 1 {M {n ■ m) - m)'^ 
8nJ (|rf + e3)^ 

1 f dV 



• I 3 I j^n^z ) +"^z 



-Rmax TT 

2 ■ sin 4> cos 4>drd(j) + 



4(47r)^ J J (r3 + e3) 





4(47r) 



(r3 + e^y 



• sin (f> dr d4>, 



where we used spherical coordinates. It follows 



^''^'^^i^^-nd^'- ^^^^ 



A2 Quadrupole Energy 

We constitute the quadrupole as a superposition of 2 dipoles at a 
distance d, with dipole moments fhi = m-zS and m2 = (— s), 
"pointing" into opposite directions. Introducing a softening length 
e, the field is given by 

B = bJip + bDip 

As the field energy does not depend on the direction of the dipole 



We define 



n± 



r± = \r±\, 
f± 
r± 

n± = n± • ei, 
and the field energy becomes 



£?B = 

where 
a = 
and 



2mi 



-2mi 



•7, 
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/■ 9 (n+n-) nlnt - 3{nl f - 3{nt f + 1 
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